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We give an explicit prescription for experimentally determining the evolution operators which 
completely describe the dynamics of a quantum mechanical black box - an arbitrary open quantum 
system. We show necessary and sufficient conditions for this to be possible, and illustrate the 
general theory by considering specifically one and two quantum bit systems. These procedures may 
be useful in the comparative evaluation of experimental quantum measurement, communication, 
and computation systems. 
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I. INTRODUCTION 

Consider a black box with an input and an output. 
Given that the transfer function is linear, if the dynamics 
of the box are described by classical physics, well known 
recipes exist to completely determine the response func- 
tion of the system. Now consider a quantum-mechanical 
black box whose input may be an arbitrary quantum 
state (in a finite dimensional Hilbert space), with inter- 
nal dynamics and an output state (of same dimension as 
the input) determined by quantum physics. The box may 
even be connected to an external reservoir, or have other 
inputs and outputs which we wish to ignore. Can we 
determine the quantum transfer function of the system? 

The answer is yes. Simply stated, the most arbitrary 
transfer function of a quantum black box is to map one 
density matrix into another, pin^Pout, and this is de- 
termined by a linear mapping £ which we shall give a 
prescription for obtaining. The interesting observation 
is that this black box may be an attempt to realize a 
useful quantum device. For example, it may be a quan- 
tum cryptography channel Jl|,H (which might include an 
eavesdropper!), a quantum computer in which decoher- 
ence occurs, limiting its performance or just an im- 
perfect quantum logic gate |5|,^| , whose performance you 
wish to characterize to determine its usefulness. 

How many parameters are necessary to describe a 
quantum black box acting on an input with a state space 
of TV dimensions? And how may these parameters be 
experimentally determined? Furthermore, how is the re- 
sulting description of £ useful as a performance charac- 
terization? 

We consider these questions in this paper. After sum- 
marizing the relevant mathematical formalism, we prove 



that £ may be determined completely by a matrix of 
complex numbers Xi an d provide an accessible experi- 
mental prescription for obtaining X- We then give explicit 
constructions for the cases of one and two quantum bits 
(qubits), and then conclude by describing related perfor- 
mance estimation quantities derivable from \- 



II. STATE CHANGE THEORY 

A general way to describe the state change experi- 
enced by a quantum system is by using quantum opera- 
tions, sometimes also known as super scattering operators 
or completely positive maps. This formalism is described 
in detail in |Q , and is given a brief but informative re- 
view in the appendix to ||. A quantum operation is a 
linear map £ which completely describes the dynamics of 
a quantum system, 



tv(£(p)) 



(2.1) 



A particularly useful description of quantum operations 
for theoretical applications is the so-called operator-sum 
representation: 



(2.2) 



The Ai are operators acting on the system alone, yet 
they completely describe the state changes of the sys- 
tem, including any possible unitary operation (quantum 
logic gate), projection (generalized measurement), or en- 
vironmental effect (decoherence). In the case of a "non- 
selective" quantum evolution, such as arises from uncon- 
trolled interactions with an environment (as in the deco- 
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herence of quantum computers), the Ai operators satisfy 
an additional completeness relation, 



E^w. 



(2.3) 



for some set of complex numbers dim- Eq.(2.2) may thus 
be rewritten as 



E^ 



(3.2) 



This relation ensures that the trace factor tr(£(p)) is al- 
ways equal to one, and thus the state change experienced 
by the system can be written 



(2.4) 



Such quantum operations are in a one to one corre- 
spondence with the set of transformations arising from 
the joint unitary evolution of the quantum system and 
an initially uncorrelated environment |?]]. In other words, 
the quantum operations formalism also describes the 
master equation and quantum Langevin pictures widely 
used in quantum optics ||,[l0|], where the system's state 
change arises from an interaction Hamiltonian between 
the system and its environment [ 111) . 

Our goal will be to describe the state change process 
by determining the operators A4 which describe £, (and 
until Section VI we shall limit ourselves to those which 
satisfy Eq.( |2.3|) ). Once these operators have been deter- 
mined many other quantities of great interest, such as 
the fidelity, entanglement fidelity and quantum channel 
capacity can be determined. Typically, the Ai operators 
are derived from a theoretical model of the system and 
its environment; for example, they are closely related to 
the Lindblad operators. However, what we propose here 
is different: to determine systematically from experiment 
what the A4 operators are for a specific quantum black 
box. 



III. GENERAL EXPERIMENTAL PROCEDURE 

The experimental procedure may be outlined as fol- 
lows. Suppose the state space of the system has N di- 
mensions; for example, N = 2 for a single qubit. N 2 
pure quantum states \tpi}{i^i\, . . . , \ipN 2 ){4>N 2 \ are exper- 
imentally prepared, and the output state £(\ipj)(ipj\) is 
measured for each input. This may be done, for example, 
by using quantum state tomography |l^-|lj]. In princi- 
ple, the quantum operation £ can now be determined by 
a linear extension of £ to all states. We prove this below. 

The goal is to determine the unknown operators Ai in 
Eq. fl2.2p . However, experimental results involve numbers 
(not operators, which are a theoretical concept). To re- 
late the Ai to measurable parameters, it is convenient to 
consider an equivalent description of £ using a fixed set of 
operators Ai , which form a basis for the set of operators 
on the state space, so that 



At 



E< 



A, 



(3.1) 



where Xmn = J2i a imO-\ n is a "classical" error correlation 
matrix which is positive Hermitian by definition. This 
shows that £ can be completely described by a complex 
number matrix, \, once the set of operators A^ has been 
fixed. In general, x w iU contain TV 4 — N 2 independent 
parameters, because a general linear map of N by N ma- 
trices to iV by iV matrices is described by N 4 independent 
parameters, but there are N 2 additional constraints due 
to the fact that the trace of p remains one. We will show 
how to determine \ experimentally, and then show how 
an operator sum representation of the form Eq.(2.2) can 
be recovered once the % matrix is known. 

Let pj, 1 < j < N 2 be a set of linearly independent 
basis elements for the space of NxN matrices. A con- 
venient choice is the set of projectors |n)(m|. Experi- 
mentally, the output state £(\n)(m\) may be obtained 
by preparing the input states \n), \m), \n + ) = (|n) + 
|m))/v / 2, and |n_) = (|n) +i\m))/y2 and forming linear 
combinations of £(|n)(n|), £(\m)(m\), £(|n+)(n+|), and 
£(|n_)(n_|). Thus, it is possible to determine £(pj) by 
state tomography, for each pj. 

Furthermore, each £(pj) may be expressed as a linear 
combination of the basis states, 



£ (Pj) = ^2^jkPk 

k 



(3.3) 



and since £{pj) is known, Xjk can thus be determined. 
To proceed, we may write 



A mPj A\=Y,^TPk: 



(3.4) 



where (3p, n are complex numbers which can be deter- 
mined by standard algorithms given the A m operators 
and the pj operators. Combining the last two expres- 
sions we have 



EEx^>b = E A 



jkPk ■ 



(3.5) 



k mn k 

From independence of the pk it follows that for each k, 



E ftjk n Xn 



A 



jk 



(3.6) 



This relation is a necessary and sufficient condition for 
the matrix \ to give the correct quantum operation £. 
One may think of x an d A as vectors, and (3 as a N 4 x TV 4 
matrix with columns indexed by mn, and rows by ij. To 
show how x may be obtained, let k be the generalized 
inverse for the matrix j3, satisfying the relation 
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X 



mn 



(3.7) 



st.xy 



Most computer algebra packages are capable of finding 
such generalized inverses. In appendix ^ it is shown that 
X defined by 



jk 



(3.8) 



satisfies the relation (3.6). The proof is somewhat sub- 
tle, but it is not relevant to the application of the present 
algorithm. 

Having determined \ one immediately obtains the op- 
erator sum representation for £ in the following manner. 
Let the unitary matrix ft diagonalize \, 



Xr, 



xy 



From this it can easily be verified that 



(3.9) 



(3.10) 



gives an operator-sum representation for the quantum 
operation £. Our algorithm may thus be summarized as 
follows: A is experimentally measured, and given (3, de- 
termined by a choice of A, we find the desired parameters 
X which completely describe £ . 



IV. ONE AND TWO QUBITS 

The above general method may be illustrated by the 
specific case of a black box operation on a single quantum 
bit (qubit). A convenient choice for the fixed operators 
At is 



Ao = I 

A~x = <?x 

A 2 = —U7y 

A 3 = cr z , 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



where the o~i are the Pauli matrices. There are 12 param- 
eters, specified by x, which determine an arbitrary single 
qubit black box operation £\ three of these describe ar- 
bitrary unitary transforms exp(i^ fc r^ffk) on the qubit, 
and nine parameters describe possible correlations estab- 
lished with the environment E via exp(z yV fe 7jfc<7j 0(7^). 
Two combinations of the nine parameters describe phys- 
ical processes analogous to the 7\ and T 2 spin-spin and 
spin-lattice relaxation rates familiar to us from classical 
magnetic spin systems. However, the dephasing and en- 
ergy loss rates determined by x do not simply describe 
ensemble behavior; rather, x describes the dynamics of a 
single quantum system. Thus, the decoherence of a single 



qubit must be described by more than just two parame- 
ters. Twelve are needed in general. 

These 12 parameters may be measured using four sets 
of experiments. As a specific example, suppose the in- 
put states |0), |1), |+> = (|0> + and |-) = 
(|0) +i \ l))/y/2 are prepared, and the four matrices 



p' 1 = £(\0)(0\) 
^ = £(|1><1|) 

P3 = £(l+X+I) 



-iSQ-)(- 

+^(l-><- 



(4.5) 
(4.6) 

-(l-i)(p[+p> 4 )/2 (4.7) 
-(l + i)(p' 1 + p' 4 )/2 (4.8) 



are determined using state tomography. These corre- 
spond to p'j = £{pj), where 



Pi = 



1 




(4.9) 



P2 = Piffx, P3 = ffxPi, and pi = a x picrx- From Eq.(3.4) 



and Eqs. (4.1-4.4) we may determine /3, and similarly p'j 
determines A. However, due to the particular choice of 
basis, and the Pauli matrix representation of Ai, we may 
express the f3 matrix as the Kronecker product (3 = A® A, 
where 



A 



I 



so that x may be expressed conveniently as 



X 



A 



Pi P2 
P3 Pi 



A. 



(4.10) 



(4.11) 



in terms of block matrices. 

Likewise, it turns out that the parameters xi describ- 
ing the black box operations on two qubits can be ex- 
pressed as 



X2 



A 2 7j'A 2 



(4.12) 



where A 2 = A® A, and p' is a matrix of sixteen measured 
density matrices, 



P 1 



Pll Pl2 Pl3 Pu 

Pix P22 P23 P24 

P31 P32 P33 P34 

P41 P42 P43 Pi4 



P. 



(4.13) 



where p' nm = £(p nm ), P™ = T n |00)(00|T m , T x = I ® J, 
T 2 = I ® ff x , T 3 = a x <g> I, T 4 = a x ® a x , and P = 
I ® [[poo + P12 + P21 + P33) ® /] is a permutation matrix. 
Similar results hold for k > 2 qubits. Note that in gen- 
eral, a quantum black box acting on k qubits is described 
by 16 fc — 4 fe independent parameters. 

There is a particularly elegant geometric view of quan- 
tum operations for a single qubit. This is based on the 
Bloch vector, A, which is defined by 



3 



p = 



I + X-a 



(4.14) 



satisfying |A| < 1. The map Eq^^l) is equivalent to a 
map of the form 



A A' = MX 



(4.15) 



where M is a 3x3 matrix, and c is a constant vector. 
This is an affine map, mapping the Bloch sphere into 
itself. If the Ai operators are written in the form 



fc=i 

then it is not difficult to check that 



(4.16) 



a lj a *k+ a *j a lk + 



N - E P a iP a ip) s jk+ 



i jp 



(4.17) 
(4.18) 



where we have made use of Eq.(2.3) to simplify the ex- 
pression for c. 



The meaning of the affine map Eq.(4.15) is made 
clearer by considering the polar decomposition |l5j ] of the 
matrix M . Any real matrix M can always be written in 
the form 



M = OS , 



(4.19) 



where O is a real orthogonal matrix with determinant 1, 
representing a proper rotation, and S is a real symmet- 
ric matrix. Viewed this way, the map Eq.( 4.15 ) is just a 
deformation of the Bloch sphere along principal axes de- 
termined by S, followed by a proper rotation due to O, 
followed by a displacement due to c. Various well-known 
decoherence measures can be identified from M and c; 
for example, T\ and T<z are related to the magnitude of 
c and the norm of M. Other measures are described in 
the following section. 



V. RELATED QUANTITIES 

We have described how to determine an unknown 
quantum operation £ by systematically exploring the re- 
sponse to a complete set of states in the system's Hilbert 
space. Once the operators Ai have been determined, 
many other interesting quantities can be evaluated. A 
quantity of particular importance is the entanglement fi- 
delity This quantity can be used to measure how 
closely the dynamics of the quantum system under con- 
sideration approximates that of some ideal quantum sys- 
tem. 



Suppose the target quantum operation is a unitary 
quantum operation, U(p) = U pU^ , and the actual quan- 
tum operation implemented experimentally is £. The 
entanglement fidelity can be defined as |l6| 

F e (p,U,E) = Y.H U]A ^\ 2 t 5 - 1 ) 

i 

= 2 Xm»tr(ff+i m p)tr(pit u) . (5.2) 

inn 

The sec ond expression follows from the first by using 
Eq.( |3.l|) , and shows that errors in the experimental de- 
termination of £ (resulting from errors in preparation 
and measurement) propagate linearly to errors in the es- 
timation of entanglement fidelity. The minimum value of 
F e over all possible states p is a single parameter which 
describes how well the experimental system implements 
the desired quantum logic gate. 

One may also be interested in the minimum fidelity of 
the gate operation. This is given by the expression, 



F = min(V>|[/ t £(|V')W'|)^7|V , ), 



(5.3) 



where the minimum is over all pure states, \ip). As for 
the entanglement fidelity, we may show that this quantity 
can be determined robustly, because of its linear depen- 
dence on the experimental errors. 

Another quantity of interest is the quantum channel 
capacity, defined by Lloyd |T^,|l8j as a measure of the 
amount of quantum information that can be sent using 
a quantum communication channel, such as an optical 
fiber. In terms of the parameters discussed in this paper, 



C{£) = max S(£(p)) 
p 



S e (p,£) 



(5.4) 



where S(£(p)) is the von Neumann entropy of the density 
operator £(p), S e (p,£) is the entropy exchange and 
the maximization is over all density operators p which 
may be used as input to the channel. It is a measure 
of the amount of quantum information that can be sent 
reliably using a quantum communications channel which 
is described by a quantum operation £ . 

One final observation is that our procedure can in prin- 
ciple be used to determine the form of the Lindblad oper- 
ator, C, used in Markovian master equations of the form 



p = £(p), 



(5.5) 



where for convenience time is measured in dimension- 
less units, to make C dimensionless. This result follows 
from the fact that Lindblad operators C are just the loga- 
rithms of quantum operations; that is, exp(£) is a quan- 
tum operation for any Lindblad operator, C, and log£ is 
a Lindblad operator for any quantum operation £. This 
observation may be used in the future to experimentally 
determine the form of the Lindblad operator for systems, 
but will not be explored further here. 
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VI. QUANTUM MEASUREMENTS 

Quantum operations can also be used to describe mea- 
surements. For each measurement outcome, i, there is 
associated a quantum operation, Ei. The corresponding 
state change is given by 



Up) 

tr(£i(p)) 



(6.1) 



where the probability of the measurement outcome oc- 
curring is j>i — tr (£,•(/?)). Note that this mapping may be 
nonlinear, because of this renormalization factor. 

Despite the possible nonlinearity, the procedure we 
have described may be adapted to evaluate the quan- 
tum operations describing a measurement. To determine 
£i we proceed exactly as before, except now we must per- 
form the measurement a large enough number of times 
that the probability pi can be reliably estimated, for ex- 
ample by using the frequency of occurrence of outcome 
i. Next, p'j is determined using tomography, allowing us 
to obtain 



£i(pj) = tv{£i{ Pj ))p' j: 



(6.2) 



for each input pj which we prepare, since each term on 
the right hand side is known. Now we proceed exactly 
as before to evaluate the quantum operation Si. This 
procedure may be useful, for example, in evaluating the 
effectiveness of a quantum-nondemolition (QND) mea- 
surement Hi. 



VII. CONCLUSION 
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APPENDIX A: PROOF OF THE X RELATION 



The difficulty in verifying that x defined by ( 3 



satis- 



fies (3.6) is that in general x IS n °t uniquely determined 
by the last set of equations. For convenience we will 
rewrite these equations in matrix form as 



X = k\. 



(Al) 
(A2) 



From the construction that led to equation (|3.2|) we know 



there exists at least one solution to equation (Al), which 
we shall call x' ■ Thus A = fix' ■ The generalized inverse 
satisfies /3k/3 = (3. Premultiplying the definition of \ by 
j3 gives 



f3x = [3 K \ 

= m% 

= A. 



(A3) 
(A4) 
(A5) 
(A6) 



Thus x defined by (A2) satisfies the equation (Al), as 
was required to show. 



In this paper we have shown how the dynamics of 
a quantum system may be experimentally determined 
using a systematic procedure. This elementary system 
identification step |2(i[ ] opens the way for robust exper- 
imental determination of a wide variety of interesting 
quantities. Amongst those that may be of particular 
interest are the quantum channel capacity, the fidelity, 
and the entanglement fidelity. We expect these results 
to be of great use in the experimental study of quantum 
computation, quantum error correction, quantum cryp- 
tography, quantum coding and quantum teleportation. 
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